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I. INTRODUCTION 

Both astronomical and cosmological data seem to require the presence of yet directly 
undetected dark matter and dark energy in the universe. The necessity for these mysteri- 
ous components occurs at distances where the gravitational interaction is not understood 
sufficiently This suspicious coincidence inspires a search for modifications of the general 
relativity at large distances. It is important to study massive and multigraviton theory 
for understanding cosmology and unification. In the linear-field theory, gravitons have the 
Fierz-Pauli (FP) type masses [1]. But there is an ambiguity in its nonlinear generalization. 
We studied thus far the linear multigraviton theory on a circle corresponding to S* 1 compact- 
ification of the KK theory with dimensional deconstruction [2]. This model is an extended 
version of Hamamoto's model [3] for a massive graviton. 

In this paper, we construct the FP Lagrangian for multigravitons associated with a gen- 
eral graph and investigate what modes of particles are included. Furthermore, we extend it 
to nonlinear theory based on the vierbein formalism [4, 5]. Nonlinear extensions of multi- 
graviton theory have been studied many authors [6]. In the present paper we focus on the 
semiclassical sector of the theory which governs the evolution of the universe; in other words, 
we will not consider nonlocal contributions and terms with higher derivatives in the possible 
complete theory here. 

The features of our model are following: (i) Gravitons as the fluctuation from Minkowski 
space-time have the FP type masses [1]. (ii) This model is based on a generalized dimensional 
deconstruction method. So, the mass spectrum in the model can be tuned more easily than 
in the KK theory, (iii) The mass term has a reflection symmetry assigned at each vertex 
and an exchange symmetry assigned at each edge of a graph. 

In this paper, beginning with graph theoretical description, we introduce the dimensional 
deconstruction [7, 8] and description of the linear theory of multigravity as the basis of our 
model in Sec. II. A nonlinear extension of the model is proposed in Sec. III. In Sec. IV, 
we consider the vacuum cosmological solutions of the case associated with the four-site 
star graph and the four-site path graph. The study on the quantum cosmological model is 
exhibited in Sec. V. Finally, we summarize our work and give remarks about the outlook in 
Sec. VI. 
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II. MULTIGRAVITON THEORY ON A GENERAL GRAPH 

A. FP on a graph 

We consider the matrix representation of the graph theory. 1 A graph G is a pair of V and 
E, where V is a set of vertices (sites) while E is a set of edges (links). An edge connects two 
vertices; two vertices located at the ends of an edge e are denoted as o(e) and t(e). Then, 
we introduce two matrices, an incidence matrix and a graph Laplacian, associated with a 
specific graph. The incidence matrix E represents the condition of connection or structure 
of a graph, and the graph Laplacian A can be obtained by EE T , where E T is the transposed 
matrix of E. By use of these matrices, a quadratic form of vectors a T Aa(= a T EE T a) can 
be written as a sum of (a t ( e ) — a D ( e )) 2 - If all dj (« = 1, 2, . . . , #V), the components of a, take 
the same value, E T a = and then Aa = 0. 

So, we consider the Lagrangian for massive gravitons h v on each vertex with the 
Stiickelberg vector fields on each edge and a scalar field (p v on each vertex: 

t 
2 



'"' [h v ^ u (EE T h^) v - h v {EE T h) v ] 



-2 [m(EA,Y + d^ v ] (d u h^ - Ph v ) ~lYl - d » A lf > ( 21 ) 



2 

veV e£E 



where Lq is the linearized Einstein-Hilbert Lagrangian: 



E 



-\d x hl v d x h v ^ v + d x h v \d u h v ^ - d„h v ^d v h v + \d x h v d x h v 



and h v = ^h^. 

This action is invariant under the following transformations: 



(2.2) 



K„ ^ h; u + d„C + d„^, A^A^ + miE'^y-d.C, ( f ) ^^ + m(ECr, (2.3) 

where £ v and ( e are parameters on each vertex and each edge respectively. The massive 
modes of vector and scalar fields are absorbed by the massive modes of graviton fields due 
to the symmetry a la Stiickelberg. 



1 Please see [9] for a brief review of application of graph theory to field theory, and textbooks [10, 11] for 
algebraic graph theory. 
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Now we examine the gauge fixing of the Lagrangian. Suppose the following gauge fixing 
terms: 

2 



L gf = - J2 \d u h v » u - -d»h v - m{EA») v - d^ v 
L 2 

vev 

- - j(E T hy + m{E T <py] 2 , (2.4) 



e&E 

then, the gauge-fixed Lagrangian becomes 

L m + L gf = X -H^[d 2 - m 2 EE T ) (h^ - iff 77^ 

+A^(d 2 - m 2 E T E)A fl + 3(p(d 2 - m 2 EE T )(f) , (2.5) 

where H^ u = h^ u + (pi]^ u . Here the indices v and e, and the notion of sum over them are 
omitted. 

In the next section, we will see that the mass spectra of fields in the Lagrangian for 
specific graphs with large number of vertices are similar to those of a five- dimensional model 
with a compactified extra space. 



B. Dimensional deconstruction 



It is assumed that we put fields on vertices or edges. An idea that there are four dimen- 
sional fields on the sites (vertices) and links (edges), dubbed as dimensional deconstruction, 
is introduced by Arkani-Hamed et al. [7, 8]. In this scheme, the square of mass matrix is 
proportional to the Laplacian of the associated graph. 




FIG. 1. The cycle graph C 60 . 

In the case of a cycle graph (a 'closed circuit') with N sites (denoted as Cn, and C 60 
is shown in Fig. 1 for example), when iV becomes large, the model on the graph coincides 
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with the five- dimensional theory with S 1 (circle) compactification. In other words, the mass 
scale of the model / over N corresponds to the inverse of the circumference of the circle: 



M 2 t = 4/ 2 (sin7r£/iV) 2 ->• M 2 = (2n£/L) 2 , (f/N -> l/L) . 



(2.6) 



The mass spectrum is given by the eigenvalues of the graph Laplacian of Cn, which can 
be expressed as 

/ 2 -1 ••• -l\ 
-1 2 -1 ••• 
0-1 2 ••• 



A = 



\ -1 



(2.7) 



For a cycle graph, the linear graviton model presented in the previous subsection coincides 
with the model proposed in Ref . [2] . The model is a most general linear multigraviton theory 
on a generic graph. 



C. Particle content in the multigraviton theory on a graph 

For this model, we investigate what modes of particles are contained. Although any graph 
is available for the model, here we consider two types, a cycle graph and a path graph 
Pn- The path graph has a simple structure like a chain, and has two ends (i — 1 and N) 
and the i-th vertex are adjacent to (i — l)-th and (i + l)-th vertices (1 < % < N). For 
example, we show C 4 and P4 in Fig. 2. The incidence matrix for C 4 is defined as 



FIG. 2. The cycle graph C4 and the path graph P4. 
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E{C 4 ) = 



( 1 

-1 






1 
-1 



-l\ 



1 



(2.8) 



\ 0-1 1 / 



and then 



E(C 4 )E(C,f 



( 2-1 -l\ 
-1 2-1 
0-1 2-1 
\-\ 0-1 2/ 

The eigenvalues of PP T are {0, 2, 2, 4} for C 4 . 

On the other hand, the incidence matrix for P4 is given by 

/ 1 o\ 



E{C,) T E{C A ) . 



(2.9) 



E(P 4 ) 



-1 1 
-1 1 

V -1 j 



(2.10) 



Thus 



E(P 4 )E(P 4 f 



( 1 -1 
-1 2 -1 
0-1 2 
\ -1 



o\ 


-1 

1/ 



(2.11) 



and 



E(P A ) T E(P 4 ) 



( 



\ 



2-1 
-1 2 -1 
0-1 2 



\ 



(2.12) 



are different in their sizes. The eigenvalues of EE T are {0, 2 — y/2, 2, 2 + y/2} and those of 
E T E are {2 — \/2, 2,2 + v^2} for P4. For P/v, it is known that the Laplacian eigenvalues 
are 4 sin 2 |^ (A; = 0, 1, ... iV — 1). If we introduce a mass scale / and consider the large N 
limit as in (2.6), we find 4/ 2 sin 2 — >■ (^) 2 where //iV = 1/L. This spectrum corresponds 
to that of the compactification on S 1 /Z 2 , where the circumference of S* 1 is 2L. 

In the multigraviton theory associated with the cycle graph Cn (#^ = N, j^E = N), 
N — 1 massive spin-two's, a massless spin-two, N — 1 massive vectors, a massless vector, 
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N — 1 massive scalars, and a massless scalar seem to be included, as seen from the gauge- 
fixed Lagrangian (2.5). The mass spectra of different spin fields are the same, except for 
zero modes. This is due to the fact that eigenvalues of EE T and ones of E T E are the same 
except for zero eigenvalues. 

However, N — 1 massive spin two, a massless spin two, a massless vector, and a mass- 
less scalar are left physically, because massive vectors and massive scalars are absorbed by 
massive spin two fields to form massive gravitons with five degrees of freedom each. 

Similarly, in the model associated with the path graph P N (#V = N, #E = N — 1), 
N — 1 massive spin two's, a massless spin two, and a massless scalar is left physically, the 
massless vector mode is absent. 

The limits of N to infinity in the cases of Cat and Pjv realize the KK theory with S l and 
S 1 /Z 2 compactification, respectively. 

III. NONLINEAR EXTENSION OF A MULTIGRAVITON THEORY ON A 

TREE GRAPH 

Now we will consider a nonlinear extension of the linear theory. Following Nibbelink 
et al. [4, 5], we introduce a useful 'tool': 



where e is the totally antisymmetric tensor. Using this expression, we have the Einstein- 
Hilbert term replacing A and B by vierbeins and C and D by the curvature 2-form. In 
addition, because the fourth power of vierbein in the angle bracket is equal to the determi- 
nant of vierbeins ((eeee) = (e 4 ) = 24|e|), this expression means that the Einstein-Hilbert 
term and the cosmological term have the similar structure. 

We now assume that the following term is assigned for each edge of a graph: 



where e\ and e 2 are vierbeins at two ends of one edge. Note that this term has a reflection 
symmetry e -H- — e at each vertex and an exchange symmetry e\ -B- e 2 at each edge. 
In the weak field limit, i.e. e± — rj + fi, e 2 = rj + / 2 , 



(ABCD) = -e abcd e^A°BlC c p Dj, 



(3.1) 




(3.2) 



((eiex - e 2 e 2 ) 2 ) = 8 (([/i] - [/ 2 ]) 2 - [(/i - / 2 ) 2 ]) + 0(f) , 



(3.3) 
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where 77 is the Minkowski metric, and [/] = tr/ for notational simplicity. This quadratic 
term corresponds to the FP mass term. 2 

On the other hand, the Einstein-Hilbert term ||e|i? contains the kinetic terms of a gravi- 
ton in the lowest order up to the total derivative: 

l -\e\R = -±d x f^d x r + d x f\d v r - d»rd u f - \d x fd x f + 0(f) , (3.4) 

and \R contains the following terms in the first order: 

l -R = -d x d x f + d»d„r + O(f ) . (3.5) 

In the case of a tree graph (a graph with no closed circuit — the path graph P N is a 
tree graph, for example), we have the nonlinear Lagrangian of multigraviton theory without 
higher derivative and nonlocal terms, 

L m = \ exp <E> \e v \ R v + ( (e o(e )e o(e ) - e t ( e )e t{e )) 2 ^ , (3.6) 

where R v is the scalar curvature associated with e v and M 2 = 3m 2 /2. The scalar zero-mode 
field $ can be identified as 4>i — 02 = • • • = 

IV. CLASSICAL COSMOLOGY OF THE MULTIGRAVITON THEORY 

Now we consider two vacuum cosmological models, associated with a four-site star graph 
and a four-site path graph respectively. Both the star graph and the line graph are tree 
graphs. The star graph consists of one central vertex and the other vertices adjacent to the 
central one. The star graph K 13 is shown in Fig. 3. The incidence matrix for K 13 is 




FIG. 3. The star graph K h3 . 



2 It is known that the asymmetric part of / can be omitted [12]. 
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E(K 1>3 ) = 



Thus 



E(K h3 )E(K 1>3 f 



I 1 1 l\ 

-10 
0-1 

\ o o -l j 



( 3 —1 —1 —1 \ 
-110 
-10 10 

V -l o o i / 



(4.1) 



(4.2) 



and 



E(K li3 ) T E(K li3 ) = 



1 2 1 
1 1 2 



(4.3) 



One can see that the eigenvalues of EE T are {0, 1, 1,4} and those of E T E are {1, 1,4} for 
the star graph K ij3 . For K\^-i, N eigenvalues of the Laplacian are {0, 1, ... , 1, N}. The 
degeneracy of N — 2 eigenvalues (= 1) is apparently due to the symmetry of the star graph. 
In the case of the star graph, the associated Lagrangian for multigravitons is the following; 



M 2 4 



L s tar = -exp<S>^2\ei\Ri + — <( e i e i - e i^if) , 



(4.4) 



i=l i=2 

where, e± is on the center of the graph. On the other hand, the Lagrangian of the case of 
the path graph is 



M 2 



Lpath = g ex P $ ^2 \ ei \ Ri + S ^ 6i6i ~ ei + iei +^ 2 ) ' ( 4 - 5 ) 
i=i i=i 

where, e\ and are on each end of the graph. 

Now let us introduce the setting for cosmology. We assume the homogeneous universe 
with a spatially-constant scalar field $(t) and the following metric; 



g^dx"dx u = -e-^dt 2 + e~ m ' A 2 {t){dr 2 + r 2 dfl 2 ) , 



(4.6) 



where A^t) (i — 1, • • • , 4) are scale factors. Then, 



(( ei e t - e 3 e 3 f) = e " 2 *(*> (<><*(*) - e a ^)(e 2a ^ - e 2a ^) , 



(4.7) 
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where Oj(i) = lnAj(i). 

We show the results of numerical calculations for the two models on the same appropriate 
initial conditions in Fig. 4 and Fig. 5. In both cases the scalar field $ behaves similarly and 
in each case scale factors repeat the increase and the decrease. The oscillation of the scale 
factors in the path graph case include more different modes than that of the scale factors in 
the case of the star graph where the degeneracy of eigenvalues exists. 

The star graph model has more symmetries than the path graph model. Therefore a lot 
of modes in the star graph are degenerate, while there is no degeneracy in the spectrum 
of the line graph. In the path graph case, increase of the number of sites gives the more 
complicated behaviors of the scale factors. On the other hand, in the star graph case, the 
symmetries are preserved even if the number of sites increases. Therefore, the behaviors of 
scale factors are much similar to those in the four-site model, essentially. 



a_l i 1 1 a_2 .; 1 1 




FIG. 4. Numerical solutions of a's and <E> in the case of the four-site star graph. 
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20 *0 f0\» 100 



FIG. 5. Numerical solutions of a's and <I> in the case of the four-site path graph. 
V. QUANTUM COSMOLOGY OF THE MULTIGRAVITON THEORY 
A. The Wheeler-DeWitt equation 

In the previous section, we have seen the oscillatory behavior in the evolution of scale 
factors. As a qualitative analysis, we only show the characteristic solutions. In fact, os- 
cillations must be dependent on the initial conditions. What are the natural conditions? 
To study the initial state, we have to consider quantum behavior of cosmology. Note that 
quantum cosmology of multigraviton theory has never been studied yet as far as we know. 

In this section we consider a minimal model based on a graph P 2 , which is shown in Fig. 6 
This model has two gravitons, 3 or two scale factors. The Lagrangian density is given by 



FIG. 6. The path graph P 2 . 
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£ = 



1 

2^ 



exp$ [\e g \Rg 



ere 



! ) + (surface terms) , 



(5.1) 



where two graviton fields are labeled by g and /. This model in this case is very similar to 
f-g gravity [17] or bigravity [18], but our model also contains a massless scalar field. 
We take the metric ansatze as follows: 

-*(-N 2 dt 2 + A 2 (m\ 



ds) 



4), 



ds 2 = e-*{-N 2 dt 2 + B 2 dQ 2 3 ) . 
These choices are equivalent to field redefinitions g^ v = e~®gffi so that 



(5.2) 
(5.3) 



1 

2^2 



+ 



| e f)|i2f) + |ef)|i2f) 

and often quoted as the choice of the Einstein frame. Here we assume that A, B and <3> 
depend only on t, and dVt\ = dx 2 + dy 2 + dz 2 . The lapse function will be set as iV = 1 after 
the calculation (by the redefinition of t). Each metric is homogeneous, isotropic, and flat in 
the Einstein frame, in the present analysis. Then the action reads 



J Ldt = j dt ^2 {e 3a (-4a 2 + $ 2 ) + e 3/3 (-4/3 2 + <i> 2 )} 



+ NM 2 e~ 2 *(e a + e f) )(e a - e 



(5.4) 



where a = In A, and (5 = ln_B. The dot indicates the derivative with respect to t. The 
conjugate variables are 



12e 3a d 



12e 3 ^d 



Nk 2 ' p Nk, 2 ' 
thus we obtain the Hamiltonian of the universe as 



H = N 



K 



K 



7T$ 



MV 



3 ( e 3/3 + e 3/3)$ 

iVV 



e a + e^)(e a -e 



^2 



(5.5) 



. (5.6) 



6 ( e 3a + e 3/3) ' * 

From the Hamiltonian, we obtain the Wheeler-DeWitt (WDW) equation for the wave 
function of the universe ty(a,(3, $) [13]. Although there are ambiguities in the ordering, we 
adopt the simple replacement of conjugate variables by the derivatives with respect to the 
corresponding dynamical variables. 4 The WDW equation for the present model is 



K 

24 



-3a 



dot' 



+ e 



-3/3 



+ 



d 2 



Q( e 3a + e W) 5$2 



+ M 2 e- 2 *{e a + e^)(e Q - e 15 ) 2 V(a, (3, $) = . 



ft\ (f , a _J>\2 



(5.7) 



3 In this case, the eigenvalues of mass are and 2M/v / 3- 

4 Another plausible choice is adoption of the Laplacian in the minisuperspace. The qualitative behavior is 
not changed by the choice of the operator orderings. 
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Now we introduce new variables x and y. They are defined as 



a + 



x = 



a- 13 



(5.8) 



Since 

d _ 1 / d d 
da 2 \dx dy 
the WDW equation (5.7) is rewritten as 

.2 ( Q2 Q2 



K 



d__l ( d_ d_ 

d/3~2 V& 7 ~ dy 



-e\ cosh3y ^dx* • dy- 



+ -x-z) - 2sinh3y- 



d 2 \ K 2 d 2 

dxdy J 6 cosh 3y d<p 2 



+ 6AM 2 e~^e 6x cosh y sinh 2 y $>(x, y, <j>) = , 
where we have also introduced = $/2 for simplicity. 



(5.9) 



(5.10) 



B. Wave-packet solutions 

To analyze the WDW equation (5.10), we assume the wave packet ansatz. The wave 
packet in quantum cosmology was originally introduced in the references [14, 15], and is 
utilized recently for various models such as in Ref. [16]. The use of the wave packet is 
crucial for the case with no special 'initial' state served as in the case with the positively 
curved homogeneous space. 

The general form of the wave function is written by 



(5.11) 



where 



K, 2 d 2 



+ 64MV 4< ^e 63; 



cosh y sinh 2 y 



ip n (x,y,<j>) = E n (x,<f>)il) n (x,y,<i>) . (5.12) 



6 dy 2 ' cosh3|/ 
We assume that x and <fi are slowly evolving variables, while y is a rapidly changing variable. 5 
In other words, we assume d x In C n <^ 1 and d^ In C n 1. 

Further we approximate the equation if y has a small amplitude. Then 

k 2 d 2 



6 dy* 



+ 64M 2 e 6x - 4 V 



ij> n (x, y, <p) = E n (x, 4>)il) n (x, y, <fr) 



(5.13) 



5 This assumption leads to a universe with the increasing mean size, which looks like our present universe. 
Some violent evolutions can occur in the very early universe, but we do not consider the possiblity here. 
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If x and are slowly-developing variables, this is no other than the equation for a harmonic 
oscillator. The differential equation 



has the solution 



r(y)-by 2 ^(y) + c^(y) = 



lMy) = exp(-- [Vb]y 2 )}l n ([b] 1/4 y), 



(5.14) 



(5.15) 



where H n is Hermite polynomial in the definition of Mathematica and ip is normalizable if 

c = c n = 2\fb i^n + , n = integer (5.16) 
Therefore the approximation gives the solution of (5.13) which leads to 



ft 



where c n is given by (5.16) with 



E n (x,(j)) = — c„, 
6 



6 = 464M 2 e 6 *- 4 ^ . 

ft 



Now the differential equation for C n becomes 

ft 2 d 2 K 2 d 2 



6 dx 2 6 dcf) 2 
and can be approximated as 



+ E n {x,<j>) 



C n (x,<P) = 0, 



dx 2 d(j) 2 
Further rewriting variables as 




ft 



C n (x, 0) = . 



3a; - 20 



30- 2x 



leads to 



Q2 Q2 - WMj^(n+l)e^ 



V5 ' 

C n (X, Z) = 0. 



ox 2 oz 2 v ft 2 V" ' 2 

Finally, separating variables as C n (X, Z) = f k (Z){p kn (X) according to 

d 2 



dZ< 



k 



fk(Z)=0, 



tells us the solution 

f k {Z) = e~ ikZ 



(fkn(X) = K 2ik/Vl 



(5.17) 



(5.18) 



(5.19) 



(5.20) 



(5.21) 



(5.22) 



(5.23) 



(5.24) 
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where K v is the modified Bessel function of the second kind with 



16M 




n + 



(5.25) 



The wave packet can be written in the form 



V = J2Mv) / dkA n (k)K 2ik/VE 

n=0 J -°° 



V5 



-ikZ 



(5.26) 



The wave function behaves oscillatory in the region X < and exponentially damps in 
the region X > 0. This is because the exponential potential 'wall' in (5.22). The amplitude 
with respect to X has a maximum peak at X ~ independently to k. Therefore the general 
wave packet, in which A{k) is taken to be a Gaussian, has a peak at X ~ 0, because other 
peaks are destructively superposed. 

The universe with X ~ is preferred in general. Even in classical solution, oscillatory y 
leads to x ~ 2/30 can be confirmed. 



C. Comparison to the case with no oscillation 

If we assume 'classically' y ~ 0, i.e., assume a = (3, WDW equation reads 

' d 2 d 2 ' 

dx 2 d<p 2 

or 

' d 2 d 2 ' 

dX 2 ~ dZ 2 

The solution of this differential equation is: 

i> = h(x - 0) + f 2 (x + <i>)=g 1 {X-Z) + g 2 (X + Z) . (5.29) 

This shows much different behaviors from the 'correct' solution of the WDW equation. No 
typical peak can be expected. This is rather trivial, but this comparison reminds us the fact 
that there is at least zero-point oscillation in any oscillatory quantum system. 



^ = 0, 



ip = 0. 



(5.27) 



(5.28) 



VI. CONCLUSION AND OUTLOOK 



We have studied the simple and Lorentz-invariant theory of multigraviton, and have 
shown typical cosmological solutions. We focused our attention on the models associated 
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with the four-site star graph and the path graph and found that vacuum cosmological so- 
lutions with the scale factors show the repeated accelerating and decelerating expansions. 
The differences between these two models were discussed from a viewpoint about symme- 
tries. By using a simplest model, we also qualitatively showed that the oscillatory behavior 
is considered as necessary in quantum universe. We should investigate more plausible and 
applicable solutions for classical as well as quantum cosmology, including usual matter. 

To this end, we should study how the gravitons and the scalar field couples to various 
matter fields. To consider various coupled fields, incorporation of supersymmetry or su- 
pergravity is also of much interest. Permitting higher derivative terms and nonlocal terms 
in the action will bring more possibilities to the completion of nonlinearity and be worth 
studying still. 

As the future works, from the mathematical point of view, it is interesting to construct 
models with the use of generic graphs, such as weighted graphs, fractals, and so on. 
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